Multifractality at the quantum Hall transition: Beyond the parabolic paradigm 
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We present an ultra-high-precision numerical study of the spectrum of multifractal exponents A, 
characterizing anomalous scaling of wave function moments dV'P'') a-t the quantum Hall transition. 
The resuh reads A, = 2g(l - q)[hQ + bi{q - l/2f + ...], with bo = 0.1291 ± 0.0002 and 6i = 
0.0029 ± 0.0003. The central finding is that the spectrum is not exactly parabolic, fei 7^ 0. This 
rules out a class of theories of Wess-Zumino-Witten type proposed recently as possible conformal 
field theories of the quantum Hall critical point. 
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Introduction The quantum Hall effect is a famous 
macroscopic quantum phenomenon [l|, [2| whose discov- 
ery gave rise to one of the most active research areas 
in condensed matter physics of last three decades. The 
plateaus with quantized values of the Hall conductivity 
are separated by quantum Hall transitions, which repre- 
sent a celebrated example of a quantum critical point in 
a disordered electronic system (for a recent review, see 
Ref.[^. Identification of the critical field theory of the in- 
teger quantum Hall transition remains a major unsolved 
problem of condensed matter physics. 

One of the key characteristics of the quantum Hall 
transition point is the multifractality spectrum govern- 
ing fluctuations of amplitudes of critical wave functions. 
Specifically, the moments of wavefunctions scale with sys- 
tem size, L, with a set of anomalous exponents Ag, 

(|^(r)p«)/(|^(r)p>«^L-A. . (1) 

(The angular brackets denote the ensemble averaging.) 
Equivalently, one often characterizes multifractality by a 
closely related set of exponents, Tq = d{q — 1) -I- Aq, or 
by its Legendre transform /(a) ("singularity spectrum"), 
aq = Tq, f{aq) — quq — Tq. Hcrc d is the system dimen- 
sionality (while d — 2 for the present case of the quantum 
Hall transition, we find it useful to keep it as d in formulas 
below), and the prime denotes the g-derivative. 

Zirnbauer'^ and Tsvelikf5','6'| conjectured that the con- 
formal theory of the quantum Hall critical point is of 
the Wess-Zumino-Witten type. These proposals implyQ 
that the spectrum is parabolic, i.e., that 7^ defined 
according to 

A, = 7,9(1 - q) . (2) 

is in fact g-independentfs], 7^ = 7. An accurate numeri- 
cal analysis of the multifractal spectrum plays therefore 
a crucial role for identification of the critical theory. 

A high-accuracy evaluation of the multifractality spec- 
trum was carried out in our earlier work Q . For this pur- 
pose, we modelled systems of a much larger size than in 
preceding works and performed averaging over a large en- 
semble of wave functions, as well as a thorough analysis of 



finite-size effects. It was found that the spectrum is close- 
to-parabolic, A, ~ 79(1 -g) with 7 = 0.262 ±0.003, thus 
showing that, if deviations from parabolicity are present, 
they are rather small (of the order of 1%). While the 
data of Ref. Q were showing some indications for such de- 
viations, they were smaller than the numerical uncertain- 
ties. The latter originate from statistical noise (limited 
size of the data set) and from finite size effects affecting 
the scaling relation ([T|) which is used to extract A, . 

The goal of the present Letter is to determine the A^ 
spectrum with an ultra-high precision and to give an ulti- 
mate answer on the question "Is the multifractality spec- 
trum of the quantum Hall transition strictly parabolic?" 
For this purpose, we improve upon the earlier numerical 
analysis in two different ways. First, we utilize a statisti- 
cal ensemble that contains approximately ten times more 
samples than the one used before Second, we employ 
a recently discovered "reciprocity relation" , 

Ag = Ai_g, (3) 

for a better control of finite-size corrections. 

Relation ([3|) implies a symmetry of the Ag spectrum 
around the point q=l/2. Consequently, an expansion of 
7g about this point has a form 

7g/d = &0 + 61(9 - 1/2)' + b2{q - 1/2)4 

To verify or to exclude parabolicity of Ag, the prefactor 
bi of the quadratic term (and possibly those of higher or- 
der terms) in Eq. ([4]) should be determined numerically. 
This is the purpose of the present work. We will provide 
numerical evidence that the corrections to parabolicity 
do not vanish. Specifically, we obtain 6o=0.1291± 0.0002 
and 61=0.0029 ± 0.0003, the non-zero bi implying that 
the parabolcity is not exact. The corresponding value of 
aQ [position of the apex of the singularity spectrum f{a)] 
is ao = rf + 70 = 2.2596 ± 0.0004. 

Method: In order to find the critical eigenstates, we 
employ the same numerical strategy that has been de- 
veloped before 0. We determine the lattice time evo- 
lution operator U for the Chalker-Coddington network 
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I with periodic boundary conditions and 
nodes. Eight eigenstates with eigenvalues clos- 



q=0.6 



est to unity are found with a standard sparse matrix 
package 13, 14, 1^ from exact diagonahzation of U. We 
study systems with L=16, 32, 64, ... , 1024 with ~ 10^ 
samples for the smallest sizes and ^ 10^ for the largest 
ones. The statistical analysis proceeds via calculating the 
average inverse participation ratios. 
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which obey the scaling law 

= Cq{N) iV-(9-i)-A,/d^ 



(5) 



(6) 



The coefficients Cq become independent of N in the limit 
N ^ oo. 

As an alternative approach, we consider the ratio 



(InF, 



■lniV+ (Inc,)', 



(7) 

whose scaling yields the exponent ag — + d. In this 
way, the exonent aq is studied directly, i. e. without in- 
voking a numerical differentiation which can significantly 
increase the error bars. In analogy with Eq. (|4]), we can 
expand around q = 1/2, 



ao + ai(<Z-^)' + .... (8) 
are related via 



The coefficients of both expansions 

ao=&o-&i/4, ai=26i-&2/2, 

The averages entering Eqs. ([5]) and ([7]) are readily ob- 
tained numerically. It is beneficial to perform the scaling 
analysis of the ratio ([7]) in addition to that of Eq. ([S]) 
for several reasons. First, the curvature of Ag is more 
clearly seen in the q-derivative, A'^. Second, the finite- 
size corrections are different in the cases of Eqs. ([5]) and 
d?]), so that an agreement between the obtained expo- 
nents provides an additional confirmation of the validity 
of the N ^ (X) extrapolation procedure. Also, the rela- 
tion 01=261—62/2 allows one to extract the coefficient 62 
of the quartic term in Eq. ^ out of parabolic fits for 7^ 
and 7q. 

Numerical results: We begin the analysis of our nu- 
merical results by verifying the reciprocity relation ([3]). 
To this end, we consider the ratio 
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9_ = 7v('^i-9-A,)/d_^ 



Cl-q{N) 



(9) 



The reciprocity relation ^ implies that the leading pow- 
ers should cancel, so that exhibits only subleading cor- 
rections in 1/N. The log-linear plot. Fig. [H upper row, 
shows that saturates in the large A'^-limit with a very 
well defined asymptotic value. Thus, we confirm reci- 
procity for the exponent spectrum of the integer quan- 
tum Hall effect, as expected. [l3| Since the exponent 




FIG. 1: Upper row: ratio of inverse participation num- 
bers Pq and Pi-g [Eq. (O] at g=0.6, 1, 1.5 (from left to right). 
The flat asymptotics indicates validity of the reciprocity re- 
lation, A5=Ai_q. Small deviations from the constant be- 
havior at largest system sizes are due to residual statistical 
noise. Lower row: analogous plots for the logarithmic deriva- 
tive Sq = (Inr,)' defined in Eq. (fTCT)l . 



relation must hold only in the asymptotic regime, we can 
draw another conclusion which is important for the sub- 
sequent analysis: the observed saturation of provides 
evidence that our numerically accessible sample sizes are 
indeed large enough in order to be able to study the true 
asymptotics. 

Similar to r^, we consider the logarithmic derivative 
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■InN+ilnCqCi^qY, (10) 
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which also saturates well inside the numerical window, 
see Fig. [l] lower row. Thus, the true asymptotics of aq 
may be studied by means of Eq. ([7]) with available system 
sizes, too. 

Having gone through important prerequisites, we now 
turn to the analysis of the main data. In order to deter- 
mine a set of relatively small exponents, Ag, to an accu- 
racy considerably better than 1%, we have developed the 
following procedure. In each panel of Fig. [5] we plot for 
fixed q a family of curves labelled by a parameter d, 



Fq{N) = p^iV^-i+^li-e)?. 



(11) 



For the particular family member, for which Fq{N) be- 
comes independent of N in the limit of large N, we can 
conclude that S="fq/d. From such a procedure we extract 
the function 7^ without having to resort to any (multi- 
parameter) fitting procedure. Similarly, by studying yet 
another family of curves, 



Fq{N) ^Lq+ lnN[l + 2q)S] 



(12) 



we have direct access also to the function 7^, see Eq. ([5]), 
without the need for numerical differentiation. 
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FIG. 2: Family of curves F,(iV) = p^iV9-i-«(<7-i)9 foj. 
5=0.6,1.1,1.5 (top; left,center,right) and g=0.4, — 0.1, — 0.5 
(bottom; left,center,right). Each data set is labelled by a 
parameter S, which increases from a minimum to a maxium 
value (given in the upper plots) in steps of 0.0004. The value 
of S for which Fq{N) is flat determines the anomalous expo- 
nent, Aq=dSq{l — q). Such saturating data sets are marked 
with filled symbols (left: A; center: O; right: □); typical 
error in the corresponding value of 5 does not exceed 0.001. 
The change in symbols (i.e. in 5) for saturating data sets illus- 
trates a 5-dependence of 7^ and thus gives direct, unprocessed 
evidence of nonvanishing quartic terms in Ag. 



The functions 7, and 7^ representing the main re- 
sult of this paper are displayed in Fig. [3l Also shown 
is /S.q/q{\ — q) as derived from the earlier evaluation of 
the exponents [9( (the size of corresponding error bars is 
indicated by dotted lines). Figure [3] clearly shows that 
the curvature in 7^, which apparently has already left its 
trace in the earlier data, now fully reveals itself thanks 
to the reduced error bars. Even more pronounced is 
the resulting structure in the derivative 7g. It is reas- 
suring to notice that the new data confirm our previ- 
ous finding for &o but provide a much better accuracy: 
60=0.1291 ± 0.0002. Our new result for the curvature of 
7^ is 61=0.0029 ± 0.0003, which is clearly non- vanishing. 
These results are in full agreement with those obtained 
by a fit to 7, (see the caption to Fig. [3]). We thus con- 
clude that, although the curvature of 7^ is numerically 
rather small (61 is approximately 50 times smaller than 
60), it is non-zero: the multifractality spectrum Ag of the 
quantum Hall transition is not parabolic. 

Surface exponents: So far a network model with torus 
geometry (i.e. without boundary) has been considered. 
Recently, it has been shown 16] that wavefunction fluctu- 
ations near surfaces exhibit their own multifractal spec- 
trum with exponents A^, defined in full analogy to 
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FIG. 3: The exponents 7q=Aq/<j(l - q) (o) and 7g=A^/(l - 
2q) (o) obtained from Fig. [2] and analogous analysis for other 
values of q. The curvature in 7, and 7, implies that the 
multifractal spectrum A, is not parabolic. Also shown are 
results of the earlier work [9j (solid line). Dotted horizon- 
tal lines indicate earlier error bars in the regime < g < 1. 
Dashed lines represent parabolic fits with 6o=0. 1291 ±0.0002, 
6i=0.0029±0.0003 (left) and oo=0.1282±0.0001, ai=0.0063± 
0.0005 (right). Combination of these data yields a rough es- 
timate of the quartic term, &2 = —0.001 ± 0.001. 



with the average (. . .)s performed over the boundary sites 
only. In general, the surface exponents A^ are not related 
to their bulk counterparts in any simple manner. We 
parametrize the surface spectrum in analogy with the 
bulk case, Eqs. ([5]), labelling the correspond- 

ing parameters by a superscript "s" . The results for 7^ 
and 7| are shown in Fig. [4] It is seen that the non- 
parabolicity of the multifractality spectrum (difference of 
7| and 7^ from a constant) is present at the boundary as 
well and is in fact considerably more pronounced that in 
the bulk. The ratio Rq=jq/'jq has a clear g-dependence, 
with a minimum at the symmetry point q ~ 1/2, where 
it takes the value i?i/2=1.434± 0.005. It is also worth 
noticing that Rq is appreciably smaller than 2, a value 
naturally expected for critical theories expressed in terms 
of a free bosonic field. 

Conclusions: To summarize, we have studied numer- 
ically the wave function statistics at the quantum Hall 
critical point. We have verified that the reciprocity re- 
lation Is]) holds and have used it to control systematic 
errors related to the finite-size effects. In combination 
with a very large size of the statistical ensemble, this has 
allowed us to reach unprecedented accuracy in determi- 
nation of the multifractality spectrum A^ , with the error 
bars reduced by almost an order of magnitude compared 
to the earlier work. The result shown in Fig. [3] reads 
A, = 29(1 - q)[bo + 6i(<7 - 1/2)2 ^^(^ „ 1/2)4 ^ ^ ^ ]^ 
with 60 = 0.1291 ± 0.0002, 61 = 0.0029 ± 0.0003, and 
62 = —0.001 ± 0.001. The obtained spectrum shows clear 
non-parabolicity, 61 7^ 0, thus excluding a broad class of 
theories of the Wess-Zumino-Witten type as candidates 
in the conformal field theory of the quantum Hall tran- 
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FIG. 4: The surface exponents 7q=AJ/g(l — q). Data are 
presented in the form analogous to the bulk plot, Fig. [S] the 
curvature is even more pronounced for the surface exponents. 
Dashed lines indicate parabolic fits with bo=0.1855 ± 0.0005, 
6f =0.022 ± 0.002 (left) and og=0.1805 ± 0.001, af=0.048 ± 
0.003 (right). The resuhing estimate for 6^ is 6^ = -0.008 ± 
0.01. 



sition. These results are corroborated by the analysis of 
the surface muthfractality. While completing this work, 
we learnt about an independent study by Obuse et al. 
17| who focussed on the surface multifractahty spectrum 
and came to the same conclusions. 

We thank A. Furusaki, I. A. Gruzberg, A.W.W. Lud- 
wig, H. Obuse, and A. R. Subramaniam for useful dis- 
cussions and for sharing their data prior to publication. 
We are also grateful to A. Tsvehk and M.R. Zirnbauer 
for instructive discussions. This work was supported by 
the Center for Functional Nanostructures of the DFG. 



[1] K. von Khtzing et al, Phys. Rev. Lett. 45, 494 (1980). 

[2] R.E. Prange, S. M. Girvin (eds.). The Quantum Hall Ef- 
fect (Springer, New York, 1992)^ 

[3] F. Evers and A. D. Mirlin. .arXiv:0707.4378l to appear in 
Rev. Mod. Phys^^ 

[4] M. Zirnbauer, 'hep-th/9905054v2. 

[5] M.J. Bhaseen et al., Nucl. Phys. B 580, 688 (2000). 

[6] A. M. Tsvelik, Phys. Rev. B 75, 184201 (2007). 

[7] While rigorously proven for a number of other Wess- 
Zumino-Witten models, the statement of parabolicity of 
Aq remains, strictly speaking, a plausible conjecture for 
the model of Ref . Ij; see also M. Bershadsky, S. Zhukov, 
and A. Vaintrob, Nucl.Phys. B 559, 205 (1999). 

[8] More accurately, the parabolicity may only hold for q < 
qc, its termination at g = gc = •■• is related to wave 
function normalization, see the review for detail. 

[9] F. Evers, A. Mildenberger and A. D. Mirlin, Phys. Rev. 
B 64, 241303(R) (2001). 
[10] A. D. Mirlin, Y. V. Fyodorov, A. Mildenberger and F. 

Evers, Phys. Rev. Lett. 97, 046803 (2006). 
[11] J.T. Chalker and P.D. Coddington, J. Phys. C 21, 2665 
(1988). 

[12] R. Klesse and M. Metzler, Int. J. Mod. Phys. C 10, 577 
(1999). 

[13] P.R. Amestoy et al. Comput. Methods in Appl. Mech. 
Eng. 184, 501 (2000); P.R. Amestoy et al, SIAM J. Ma- 
trix Anal. Appl. 23, 15 (2001). 

[14] J. W. Demmel et al, SIAM J. Matrix Anal. Appl. 20, 
720 (1999). 

[15] R.B. Lehoucq, D. Sorensen and C. Yang, ARPACK Users 

Guide (SIAM, Philadelphia, 1998). 
[16] A. R. Subramaniam, I. A. Gruzberg, A. W. W. Ludwig, 

A. D. Mirlin, A. Mildenberger and F. Evers, Phys. Rev. 

Lett. 96, 126802 (2006). 
[17] H. Obuse, A. R. Subramaniam, F. Furusaki, 

I. A. Gruzberg, and A. W. W. Ludwig, preprint 

larXiv:0804.2409l 



[*] Also at Petersburg Nuclear Physics Institute, 188300 
St. Petersburg, Russia. 



